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We introduce a technique to decompose the scattered near field of two-dimensional arbitrary 
metaatoms into its multipole contributions. To this end we expand the scattered field upon plane 
wave illumination into cylindrical harmonics as known from Mie's theory. By relating these cylin- 
drical harmonics to the field radiated by Cartesian multipoles, the contribution of the lowest order 
electric and magnetic multipoles can be identified. Revealing these multipoles is essential for the 
design of metamaterials because they largely determine the character of light propagation. In par- 
ticular, having this information at hand it is straightforward to distinguish between effects that 
result either from the arrangement of the metaatoms or from their particular design. 
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I. INTRODUCTION 

Metamaterials may be understood as a kind of arti- 
ficial matter that allows to control the mould of light 
predominantly by the geometry of their building blocks 
rather by their intrinsic material properties. Fascination 
arose since these building blocks, commonly called the 
metaatoms, can be designed to allow for propagation ef- 
fects inaccessible in natural materials. To simply describe 
the optical action of metamaterials, effective properties 
are in most cases assigned that are retrieved from the 
optical response of an ensemble of metaatoms instead 
from these individual atoms themselves. Fo r this pur- 
pose, single layei'i' or bulk metamaterial^^El are usually 
treated as black boxes to which effective properties are 
assigned with the only purpose to reproduce scattering 
data like reflection and transmission coefficients. These 
data do not provide sufficient insights into the physics of 
metamaterials since their rational design usually aim s at 
evoking a certain multipolar scattering respons^^. Al- 
though metamaterials cover a wide range of structures 
at present, media with a magnetic responsd^ at optical 
frequencies are particularly appealing since they do not 
exist in nature. The pertinent metaatoms shall than pos- 
sess a strong magnetic dipole moment. To achieve this, 
the metaatoms are typically ring-shaped*^, resulting in a 
ring-like current distribution at resonance. Hence, the 
optical response contains a strong magnetic dipole field 
contribution leading to an appreciable dispersion in the 
effective permeability. Such understanding of metama- 
terials is very versatile as it provides the possibility to 
optimize metaatoms for different spectral domains. 

The advantage of understanding the optical response 
in terms of multipole scatt ering is furthermore proven 
by various theoretical workaSEEHSl, Jt was shown that 



the optical response of single metamaterial layers as 
well as of bulk metamaterials can be described by 
assuming induced multipole moment densities up to the 
second order. From the multipolar contributions of the 
field scattered by the metaatoms it is even possible to 
directly assign effective material parameters. Even the 
effects of disorder can be studied and understood in 
terms of a multipole analysis as shown theoretically and 
experimentally in Ref. 11' Hence, a detailed quantitative 
theoretical study of the optical response of the single 
metaatom is in order and in most cases sufficient to 
deduce the optical response of an ensemble of these 
entities. Although important, such a contribution is 
currently missing. Whereas first attempts are reported 
in literaturd^ the analysis was usually restricted to the 
far-field scattering response. Then, either by optimizing 
the magnitude of the different multipolar contributions 
to match a certain angular scattering response or by 
probing for the scattering strength in certain directions 
where some multipole moments do not radiate, the 
multipolar response can be revealed. Nevertheless, it 
remains an open question how unique the assignments 
based on the far fields are. 

In this contribution we develop a rigorous method to 
analyze the scattered near-field of individual metaatoms 
that permits disclosing their multipolar scattering con- 
tributions. The key ingredient is an expansion of the 
scattered field of the metaatoms upon plane wave illumi- 
nation into cylindrical harmonics, i.e., we are restricting 
the current analysis to two-dimensional metaatoms. By 
relating these cylindrical harmonics to the field of Carte- 
sian multipoles, it is possible to calculate their spectrally 
resolved amplitudes. With this method at hand we will 
subsequently investigate the multipole contributions to 
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the scattered field of two prominent and frequently stud- 
ied metaatoms providing artificial magnetism; namely 
the split-ring resonator (SRR) and the cut-wire pair 
(CW). It is shown that the scattering response contains 
contributions of electric and magnetic dipoles, but also of 
an electric quadrupole. It resonates simultaneously with 
the magnetic dipole and its contribution is much stronger 
for the CW when compared to the SRR. 

Although we are focussing here only on the analysis of 
two specific metaatoms, the present technique is general 
and can be applied to various other metaatoms as well. 
Perspectively it will permit to design metaatoms with 
specific predefined multipolar contributions to the scat- 
tered field and represent a tool to distinguish between 
properties emerging from the periodic arrangement of 
metaatoms or from the specific metaatom scattering re- 
sponse. This will be of particular importance for the pre- 
diction of effective properties of self-organized, bottom- 
up metamaterials which might not allow for a perfectly 
periodic metaatom arrangement. 



II. MULTIPOLE EXPANSION OF THE 
TWO-DIMENSIONAL SCATTERED FIELD 

In order to reveal the multipolar character of the field 
scattered by an arbitrary shaped metaatom, we will de- 
velop a method to expand its scattered field into multi- 
pole fields. Since we focus here on two-dimensional struc- 
tures, cylindrical harmonics are an appropriate system of 
eigenf unctions. As soon as we have the corresponding ex- 
pansion at hand, we will show the equivalence between 
these eigenf unctions and the Cartesian multipole fields 
for cylindrical sources. Hence, this expansion will allow 
for a direct calculation of the desired multipole coeffi- 
cients. 

We will start by briefly deriving the mathematical back- 
ground to decompose the scattered field of a single 
metaatom into multipole fields. As known from Mie's 
theorjE^the translational invariance in z-dircction allows 
to separate the general vectorial scattering problem into 
the two scalar cases of TE- and TM-polarization. For 
light propagating in a linear, homogenous, isotropic, lo- 
cal medium the tangential fields satisfy the scalar wave 
equatioiPS 
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with fc(w)^ = uP'jc?. In Eq. ([T]) Fz{x,y) denotes the 
tangential component of either the magnetic field for TM 
polarization (magnetic component out-of-plane) or the 
electric field for TE polarization (electric field component 
out-of-plane) . 

Casting Eq. ([ij into polar coordinates with x = 
Rcos{4>), y = i?sin(0), R = \/x^ + y'^ we obtain [F^ — 
F.{R,<j>)] 
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The general solutions to this equation are radially depen- 
dent Bessel functions with an azimuthally varying phase 
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Eq. ([3]) corresponds to the linearly independent solutions 
of Eq. ([2| which consist of combinations of Bessel func- 
tions of the first [Jm(fci?)] and the second kind [^^(fci?)] 
denoted as Z„i multiplied by an angular function. Finally 
the scattered field from any cylindrical object centered at 
the origin and subject to the Sommerfeld radiation con- 
dition at infinity reads as 
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where H^\kR) = J™(/ci?) + iY„,{kR) are Hankel func- 
tions of the first kind and are the expansion coeffi- 
cients, termed Mie scattering coefficients of the respective 
expansion order m. Note that any field outside a virtual 
cylinder that entirely contains the scattering object can 
be expressed by Eq. Q since these eigensolutions form 
an orthogonal and complete set of eigenfunctions. By ex- 
panding the field scattered by an arbitrary particle in this 
base we can determine the contribution of the respective 
expansion order m to the total scattered field. 

The individual scattering order m always consists of 
two parts, namely the contributions from m and —to in 
Eq. Q. Hence the scattered field (Fz^s) summarizing 
these orders can be rewritten 
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where we used [H^^ = H'l^]. Now Eq. (|5|) can be 
considered as an usual series expansion with respect to 
eigenfunctions of increasing order. These eigenfunctions 
-0™ {R, (f>, k) are simply products of a radially and an 
azimuthally varying function. Regarding the azimuthal 
terms it can be seen from Eq. ([5) that they split into 
two contributions that are Tr/2m phase-shifted indicated 
by the superscripts ±, respectively. This is related 
to the fact that each expansion order is composed 
of two physically identical but azimuthally rotated 
contributions of order to. For realistic metaatoms the 
fundamental azimuthal part can be selected regarding 
the symmetries of the nanostructure; in general both 
linear independent contributions have to be considered 
for the respective order to. 

In order to decompose the scattered field of an arbi- 
trary metaatom into the derived set of eigenfunctions 
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FIG. 1: (Color online) The modulus of the magnetic field distribution (TM polarization) described by the eigenfunctions 
according to Eq. ([5| for the first three expansion orders. We have plotted the eigenfunctions accounting for the cos{m(j)) 
angular distributions only. The results for the sm{m(j>) terms would result in a rotation by an angle of n/2m. By comparison 
with the field distributions of point multipoles, which are shown in the gray-scaled insets, it can already be anticipated that 
the eigenfunctions can be identified with multipole fields. Here m = (a) corresponds to a magnetic dipole, m — 1 (b) is 
related to the electric dipole, while the electric quadrupole is associated with m = 2 (c). 



[Eqs.([5])] the Mie coefficients a™ are required. As usual 
they are obtained by evaluating the overlap integral 
between the eigenfunction ^p^{R, cj), k) and the respective 
field component Fz {R, 0) of the individual scattering 
object 



With Eq. (|6| it is possible to rigorously determine 
based on the field overlap calculated on an annulus with 
the two radii i?2 and i?i . Here the orthogonality of the 
eigenfunctions ip^{R,(t),k) was exploited, i.e. replacing 
Fz{R,(j)) by any eigenfunction of the order I such that 
FziR,(b) = MR,<l^,k) yields 
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By exploiting the orthogonality of the azimuthally vary- 
ing part of ipi^{R, (j), k), it suffices to evaluate the over- 
lap integral for a fixed value of R rather than an an- 
nulus to obtain af^^. However, it turned out that for 
numerically (or potentially experimentally) determined 
scattered fields of realistic metaatoms, the annulus inte- 
gration is more stable and it was hence retained. This is 
not a numerical inaccuracy but it is rather attributed to 
the discrete mesh where the numerical data is available. 
This mesh is not aligned with a cylinder surrounding the 
object and an appropriate interpolation is required. Al- 
though a sufficiently fine grid improves the stability this 
issue can be circumvented if the annulus integration is 
performed. 

Although, it suffices to know the amplitudes a^j to 
fully describe the scattering response, they do not pro- 
vide the physical insight as provided by Cartesian mul- 
tipole moments. Therefore, they need to be unambigu- 
ously related. Since we have restricted ourselves to a 



two-dimensional configuration, we have to relate the scat- 
tered field of a line source exhibiting any of the relevant 
electric or magnetic multipole moments to the respec- 
tive cylindrical eigenfunctions. Starting with the well- 
known expressions for clcctrodynamic point multipoles 
in Cartesian coordinateS^^*^ we derive below how they 
are related to ipm{R,4'ik). This is performed for the 
three lowest orders by integrating the point multipoles 
in Cartesian coordinates along the z-direction to disclose 
the radiation pattern of a line source exhibiting such 
point multipoles. 



The vector potential of monochromatic electromag- 
netic fields originating from localized oscillating currents 
for the two lowest orders of a multipole expansion reads 
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FIG. 2: (Color online) The illumination conditions, the orientation, and the definition of the geometrical parameters of the CW 
(a) and the SRR metaatom (b). The calculated Mie coefficients, related to the magnetic dipole (MD), the electric dipole (ED) 
and the electric quadrupole (EQ) for the CW (c) and the SRR (d) metaatom. Finally, the scattered magnetic fields for the 
magnetic (e) and the electric resonance (f) of the CW and the SRR (g,h) are shown, respectively. The gray-scaled insets show 
the exact multipole magnetic field distributions as in Fig. [l]to underline the similarities to the exact scattering field patterns. 
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The first term in Eq. ([s]) accounts for the electric dipole 
moment where r is the length of the three dimensional 
radius vector (r — x'^ + + z^) and n = r/r is the 
normal vector. The second term consists of the electric 
quadrupole and the magnetic dipole contributions, both 
representing second-order moments in the multipole ex- 
pansion. After some algebra the electric dipole term can 
be written 
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with p(r) = J cfir'r' p{r') being the electric diploe mo- 
ment. For electric dipole moments localized in the {x,y) 
plane the z component of the magnetic field can be cal- 
culated as 
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Since we are interested in the two-dimensional represen- 



tation of the fields we integrate Eq. ( 10 ) along the z axis 
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As it can be easily verified, this electric dipole field co- 
incides with the field in Eq. ^ ior m — 1. Hence, the 
eigenfunction of the order m = 1 represents the elec- 
tric dipole contribution. Likewise one can show that the 
magnetic dipole and the electric quadrupole contribution 
correspond to m = and m = 2, respectively. A de- 
tailed derivation for the multipoles can be found in the 
appendix. Based on these results we can conclude that 
simply by calculating the Mie coefficients we have the 
multipole coefficients of the scattered field at hand. 

In Fig. [l] the magnetic field patterns for m — 0, 1, 2 
are shown. Obviously they correspond to the known field 
distributions for the magnetic dipole (m = 0) , the electric 
dipole (m = 1) and the electric quadrupole (m = 2) as 
is clear by comparison with the field patterns of point 
multipoles (see insets in Fig. [I]). 

Having finished the analytical treatment, we are going 
to apply the results to exemplarily reveal the multipole 
scattering contributions for two prominent metaatoms in 
the following. 
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III. MULTIPOLE SCATTERING OF 
METAATOMS 

For the application of the technique developed above, 
we selected the CW and the SRR geometry. We em- 
phasize that this is not a necessary restriction because 
arbitrary structures can be investigated too with the 
developed formalism. We selected these two structures 
since several modified, more complex metamaterials are 
composed out of these basic plasmonic entitied^li^l ^j^g 
investigated mctaatoms for both scenarios are shown in 
Fig. [2ja,b). The CW structure has a wire distance of 
h — 60 nm, a wire thickness of i = 20 nm and a width 
of & = 100 nm. To evaluate the overlap integral the 
annulus radii were set to Ri = 70 nm and i?2 = 80 nm. 
In order to observe the localized eigenmodes at similar 
spectral positions we used the same dimensions for the 
SRR with an additional connection {d — 20 nm) of both 
wires [Fig. 2lb)]. As a material for both metaatoms we 
selected gold^^ embedded in vacuum. The structures 
were illuminated according to the conditions as shown 
in Fig. [2]ja,b) with monochromatic plane waves. 

In order to calculate the electromagnetic near fields 
we applied the finite element method (FEM)^^. By cal- 
culating the scattering patterns for both structures and 
performing the field overlap calculations according to Eq. 
([s]) for the first three orders one obtains the Mie coeffi- 
cients as shown in Fig. [2f^c,d). For both metaatoms the 
observed low energy resonance peaks occur for both sec- 
ond order multipole contributions (electric quadrupole 
and magnetic dipole) at the same frequency (wavenum- 
ber), while the high frequency resonances are associated 
with the electric dipole modes. Thus, the usual argument 
that the electric dipole contribution of the two currents, 
oscillating tt out-of-phasc in the wires perpendicular to 
the propagation direction, annihilate and hence the next 
higher order multipole moments prevail, is confirmed by 
these results. Furthermore, the quantitative contribution 
of each multipole moment is clearly revealed. For the 
CW structure [Fig. ^c)] the electric quadrupole contri- 
bution is much stronger than that for the SRR. This is 
due to the different symmetry of both metaatoms. The 
shortcut of both wires essentially prevents the excitation 
of a quadrupole moment. Hence, the SRR's optical re- 
sponse is mainly governed by the electric and the mag- 
netic dipole moment. 

This argumentation is supported by considering the 
near fields at resonance. Comparing the pure multi- 
pole fields (Fig. IT| with the scattered fields of the two 
metaatoms [Fig. ^e-h)] it becomes obvious that the 
magnetic field of the high frequency resonance (electric 
resonance) compares nicely to that of an electric dipole. 
The scattered magnetic fields for the low frequency reso- 
nance (magnetic resonance) for the CW structure shows a 
combination of centro-symmetric magnetic dipole fields 
in between the wires and the electric quadrupole fields 
around. This is indicated by the fourfold patterns in the 
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FIG. 3; (Color online) (a) The SRR geometry together with 
the three positions representing the selected origins for the 
multipole expansion. "A" is associated with the center of 
symmetry of the corresponding CW structure, while "B" 
and " C" represent additional positions out of the metaatoms 
center. The calculated Mie coefficients for the magnetic 
dipole (MD)(b),the electric dipole (ED) (c) and the electric 
quadrupole (EQ) (d). For completeness both electric dipole 
moments along the x (a^) and the y axis (a]~) are shown. 



iso-surface lines in Fig. [2](e) outside the CW geometry. 
Differing from these patterns the SRR exhibits near field 
features that manly attributed to the magnetic dipole 
radiation pattern [Fig. [2]jg)]. In addition a weak back- 
ground contribution of an electric dipole tilted by 45° can 
be observed in Fig. [2|^g), caused by the superposition of 
weak SRR electric dipole moments in x and y directions 
[e.g. see Fig. 

Note that with the knowledge of the scattering response 
of the isolated metaatoms other physical observables be- 
come easily accessible. Most notably various other im- 
portant metamaterial properties, i.e. the effective cross 
sections^^ or the isolated polarizabilitie^^ that can be 
obtained by this Mie theory based formalism straightfor- 
wardly, can be analytically calculatecP^. 

Finally, the dependence of the multipole coefficients 
on the choice of the origin is investigated. From elec- 
trostatics it is known that the first non-vanishing multi- 
pole moment does not depend on the choice of the origin. 
For optical fields this is not valid anymore. However, for 
the scattering response of very small nanoparticles where 
the quasi-static approximation is valid, this dependency 
is expected to be negligible for reasonably chosen ori- 
gins. Exemplarily the multipole expansion is performed 
for the scattered field of a SRR with smaller dimensions 
(h = 70 nm, 6 = 70 nm, i = 25 nm, d = 20 nm) to prove 
the quasi-statically expected behavior for three different 
origins. The corresponding results are shown in Fig. [3] 
Clearly the electric dipole moments are almost constant 
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for different origins, whereas the magnetic dipole as well 
as the electric quadrupole moments may appreciably de- 
viate, in particular for the second resonance. Of course, 
the multipole moments will strongly deviate from the re- 
sults shown here if the origin is placed far outside the 
metaatoms, but then any multipole expansion becomes 
meaningless, too. Physically, the origin should be chosen 
such that all higher order multipoles beyond the second 
order are strongly suppressed. 

Nevertheless from these results we can conclude that 
for reasonably chosen origins, i.e. placed close to the 
center of mass as suggested in Ref. [10] the multipolar 
character of the scattered field can be revealed consis- 
tently. For calculating effective material parameters on 
the other hand this origin dependence has to be kept in 
mind^^. 



IV. SUMMARY 

In summary, we presented a novel technique based on 
Mie theory to reveal the contribution of various multi- 
pole excitations to the near-field scattering pattern of 
metaatoms. For the sake simplicity we restricted our- 
selves to two-dimensional metaatoms. Since Mie theory 
has been originally developed for three-dimensional ob- 
jects (spheres) the approach can be also extended to- 
wards three-dimensional metaatoms. With the presented 
formalism we revealed the quantitative contributions of 
multipole moments to the plasmonic eigenmodes of the 
cut-wire and the split-ring resonator structures. As pre- 
viously phenomenologically interpreted, due to the sim- 
ilarities of the multipole moments and the electron dy- 
namics, we could rigorously confirm the excitation of up 
to second order multipoles as the dominant scattering 
contributions. We anticipate that such a multipole ex- 
pansion of the scattered field is a genuine approach to 
optimize metaatoms with regard to a predefined scat- 
tering response. It will be very helpful to identify the 
effect of periodicity or its absence, in case of metama- 
terials fabricated with self-organization techniques, with 
respect to the effective MM properties. Beyond possible 
applications in the field of metamaterials, the presented 
approach may also provide guidelines in the design of op- 
tical nanoantennas to achieve a desired radiation char- 
acteristics to be matched to an, in principle, arbitrary 
source. 
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Appendix: The connection between 3D and 2D 
radiated multipole fields 

1. Electric dipole 

We start our considerations with the vector potential 
for the first two expansion orders Eq. ([s]) 
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At first we consider the first order term on the right hand 
side of Eq. (A.l I which represents the electric dipole con- 



tribution. Applying the continuity equation and requir- 
ing that the current density vanishes at infinity, we arrive 
at the well-known vector potential of an electric dipole 
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According to Maxwell's equations we obtain the z- 
component of the magnetic field by applying the curl 
operator to the potential [Eq. (A.2)] 
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In order to obtain the associated two-dimensional fields 
we integrate the three-dimensional fields along the z axis. 
This is equivalent to the transition from point multi- 
poles to line multipoles, similar to the transition from 
the three-dimensional to the two-dimensions Green's 
function.lani 
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The integral can be carried out with the help of Eq. 8.421 
iiP 
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Now the recursion formula for Hankel's functions of the 
first kind can be applied to calculate the remaining 
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which aUows to calculate the required magnetic fields as 

B',^jR,cl>) = c.^[cos(0K-sin(0K]i/f'(fci?). 

(A.7) 

According to Eq. ([s]) the associated magnetic field for 
m = f is given by 

Fz,s = a+'tP+{R,(j),k) + a:[^p:[{R,(j),k), 

= [a+cos{(j))-a^sm{(l))]H[^\kR). (A.8) 

It becomes obvious that the coefficients represent two 
electric dipoles rotated by 7r/2, since a comparison of the 
coefficient between Eq. (A.7) and Eq. (A.8) yields 
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2. Magnetic dipole 



The second order expansion term in Eq. (A.l) can 
be split into two contributions where one is associated 
with the vector potential of the electric quadrupole and 
the other one with the magnetic dipole. Hence, both 
quantities represent second-order moments 



A(r) 
Aoq(r) 



Mo e"" / I 
47r r \r 
Amd(r) + Aeq(r) 



-ik] d-^r'j(r')(n- r') 



Mo 


^tkr /I 


47r 


r \r 




/>oo 


X ^ 






— OO 



ik 



y ^dV-[(n.r')j(r') + (n-j(r'))r'], 



Amd(r) 



Mo e_ 
47r r 



ikr 



ik 

r 



1 

U dVr'xj(r') 

^ J-oo 



Mo e'-' 
47r r 



— ik ] m X n. 



(A.IO) 



Currents in the (a;, y) plane only induce a magnetic dipole 
moment m pointing into z direction 
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Using the identity 
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Eq. ( A.II ) can be simplified yielding 
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In order to get the magnetic field Eq. (A. 3) can be ap- 
plied which has to be integrated along the z axis as for 
the electric dipole before 



B 



md, 



Mo"^2 / 
47r \dx'^ ' dy'^ 



Q2 ^^ ^tkr 



(A.14) 



Eq. 


A.14 


(A.5 





.Mom(^)™ 



2 7.2 



H[^\kR)- H^^'>{kR) 



kR 



(A.15) 

which can be simplified using the theorem 8.473, 3iiP 



H, 



(1,2) 



(z) = -^H['''\z)-Hi'^'Hz), (A.16) 



resulting in the expected radially dependent expression 
for the scattered magnetic field 



B. 



2D 

md, 



(i?,< 



• Mo 



m,k 



H^^\kR). (A.17) 



For comparison Eq. ([s]) with m = yields 
F,,, = aoHl^\kR). 



(A.I8) 



ComparingEq. (A.I7)andEq. ( A. 18 ) we see that the co- 
efficient flo again is directly proportional to the magnetic 
dipole moment, similar to the electric dipole moment be- 



fore [Eq. (A.9)]. Here the second term connected to 
a rotated magnetic dipole by 7r/2 vanishes, as expected 
for the radial symmetric radiation pattern of a magnetic 
dipole. 



3. Electric quadrupole 

The second term connected to the electric quadrupole 



moment in the vector potential of Eq. (A.IO) can be 
simplified using Gauss' law and the continuity equation 

tjSl 



A fr) - ^ — 



— ik 



/OO 
dVr'(n-r')p(r'). (A.I9) 
-OO 
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FIG. 4: (Color online) The exact field patterns for the z component of the magnetic field for the magnetic dipole (a), the electric 
dipole (b) , and the electric quadrupole (c) for carrier dynamics in the {x, y) plane. The two-dimensional field distributions 
below the three-dimensional ones show the three-dimensional fields for the respective multipole moment integrated along 
the z axis. The moduli of these two-dimensional fields correspond to the insets in Fig. [ija-c) and precisely agree with 
ipmir, (f), m £ {1, 2, 3}, respectively as shown in the derivation below. 



In the following we restrict the carrier dynamics to the 
(x, y) plane that might be expressed by 



In a final stepEq. (A. 14) can be integrated by applying 
Eq. (TOl) 



N 



p(r') = J2 lo^SW - x^m\^ ~ y^{t)W). (A.20) 



• Mo 



Qxy 



We now rewrite the integral kernel of Eq. (A. 19) 

3 3 »oo 



{Qyy QxxJ 



dxdy 



dV'r'(n . r')p(r') = E E / dVx^n^Xlpir'), 

(A.21) 



(A.25) 



Qji{n)ej. 



Here we have used Q{n) 



N I xln^ XayaTly 

Q(n) = VaXarix yluy 

a=l \ 

Qxx^x Qxy^y 

Qxy^x Qyy^y 



(A.22) 











where the Qg are the primitive sy mmetri c quadr upole 
tensor entriePl. Substituting Eq. (|A.22|), (|A.2l|) into 



potential 



Eq. (|A.19| and applying Eq. (|A.12|) we obtain the vector 

d 



Acq(r) 



d 



d 



Cxx o ' Qxy r^ 

dy 



dx 



+ey ( Qyx + Qyy -7^ 



This equation can again be u sed to calculate the respec- 
tive magnetic field [Eq. (A.3|] 



Bcq,.(r) = - 



Mo 



92 



vv 



Qx 



d 



^ikr 



oxay ayox 



which can be simplified using the recursion formula Eq. 
(A.6I for the second derivatives 



-»^{Ox,[cos2(0)-sin2(0)] 

o 

+ [Qvv - Qxx) cos((/)) sin((/.)} i/f ^ {kR). 

(A.26) 



Finally we end up with the magnetic field 



= -i— — [(5xyCos(2(/)) 
16 



{Qyy~Qxx)sm{2^)]Hi^\kR). 

(A.27) 



■ (A. 23) ^j^g evaluation of Eq. ([5| for m = 2 yields 



F,^s = - [at cos(20) + ia^ sin(2(/))] H^/\kR) 



(A.28) 



Again the coefficients 02 are proportional to the electric 
(A. 24) quadrupole moment entries under consideration. 
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